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Abstract 

From eigensolutions of the harmonic oscillator or Kepler-Coulomb Hamil- 
tonian we extend the functional equation for the Riemann zeta function and 
develop integral representations for the Riemann xi function that is the com- 
pleted classical zeta function. A key result provides a basis for generalizing the 
important Riemann-Siegel integral formula. 

Key words and phrases 

harmonic oscillator, Mellin transformation, Hermite polynomial, hypergeometric se- 
ries, functional equation, theta function, Riemann xi function 



PACS classification numbers 

02.30.Gp, 02.30.Uu, 02.30.-f 

AMS classification numbers 

33C05, 42C05, 44A15, 44A20 

Author contact information 

fax 303-273-3919, e-mail mcoffey@mines.edu 



1 



Introduction 

The quantum harmonic oscillator and Kepler-Coulomb problems are of endur- 
ing interest to mathematical physics and are in a sense the same problem, for there 
are many ways to transform one to the other. For example, the 4-dimensional (8- 
dimensional) harmonic oscillator may be transformed to a 3-dimensional (5-dimensional) 
Coulomb problem [T9~ll2~T]. Many more mappings are realizable, especially for the cor- 
responding radial problems [19J. Therefore, the eigensolutions are very closely related 
for these two problems with central potentials. 

We recall that the eigensolution of the fundamental quantum mechanical prob- 
lems of the harmonic oscillator and hydrogenic atoms contain Hermite or associated 
Laguerre polynomials, depending upon the coordinate system used and the spatial 
dimension (e.g., [2H HD])- Two quantum numbers appear for indexing the energy 
levels and the angular momentum. These wavefunctions have a wide variety of ap- 
plicability, including to image processing and the combinatorics of zero- dimensional 
quantum field theory [TO]. Very recently additional analytic properties of these "quan- 
tum shapelets" have been expounded [10J. The one- dimensional Coulomb problem 
has recently reappeared as a model in quantum computing with electrons on liquid 
helium films [26l [25] . In addition, given the self- reciprocal Fourier transform prop- 
erty of Hermite polynomials, there are several applications in Fourier optics jH [16]. 
Hermite and Laguerre polynomials are also important in random matrix theory, es- 
pecially for determinantal processes [HI [23]. There, the joint probability density 
function of the eigenvalues of a matrix from one of the Gaussian invariant ensembles 
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is proportional to a Vandermonde determinant and the orthogonal polynomials form 
the kernel giving the n-point correlation function. 

The harmonic oscillator Hamiltonian may be given a group theoretic interpretation 
via the Weil representation of SL 2 (R) [20] and such theory together with that of 
quantum mechanical commutation relations figures prominently in constructing theta 
functions [291 U\- From the Mellin transform of suitable theta functions, one may 
represent completed zeta functions and one of our main results demonstrates this. 
Additionally, we present representations of the Riemann xi function that offer the 
possibility to generalize the very important Riemann- Siegel formula in the theory of 
the Riemann zeta function [T2l [28] . 

The work of Bump et al. [HI El EH] on the Mellin transforms of Hermite and as- 
sociated Laguerre functions has created significant interest since the zeros of these 
functions lie only on the critical line Re s = 1/2. A complementary point of view 
is possible within the theory of special functions and in particular, when the Mellin 
transforms are written in terms of the Gauss hypergeometric function, well known 
transformation formulae yield functional equations, reciprocity laws, and other prop- 
erties [H] . As an example, we have very recently shown that the Mellin transform 
of the Laguerre function £"(x) = x a ^ 2 e~ x ^ 2 L^(x), where a > — 1 and is the asso- 
ciated Laguerre polynomial [21 [151 [22] , is given by M£(s) = 2 s+a / 2 T(s + a/2)P«(s) 
where T is the Gamma function and P"(s) = <yl+ ^ )n 2F\{—n, s+a/2; a+1; 2), with p F q 
the generalized hypergeometric function and (a) n = T(a + n)/T(a) the Pochhammer 
symbol. In fact, the polynomials P™ are closely related to the symmetric Meixner- 
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Pollaczek polynomials P( x \x,n/2) [3l [18], and the latter are very useful in random 
matrix theory |14j . 

Keating [T7] showed how to generalize Riemann's second proof of the functional 
equation of £(s) by using Mellin and Fourier transforms of Hermite polynomials. In 
the next section we illustrate the hypergeometric function point of view of those 
results. The succeeding section contains our key result. We extend earlier work 
[5] such that the representations given for the xi function could provide a method 
to generalize the Riemann-Siegel integral formula of complex analysis. Again we 
emphasize that these representations are built from the eigensolutions of self-adjoint 
operators, these operators being the distinguished harmonic oscillator or Coulomb 
Hamiltonians. 



In Ref. [T7] Keating generalized the Jacobi inversion formula, applying Poisson 
summation to Hermite functions TC n (x) = x e H n (ax)e~ cx2 , where H n is the Hermite 
polynomial, a ^ 0, c = a 2 /2, and £ = 0, 1, 2, . . .. Let the 6 function u n j(t) be given 



A family of zeta functions 



by 



oo 




(1) 



m=l 



Then we have 



Theorem. [32] For n = 2q, n + 1 = 2p, p, q, I = 0, 1, 2, . . . , Re s > 1 + £, 
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where 



and 



H q (s) 



H q (s) 



(3) 



(4) 



fc f A;!(2g-2A;)! V27 V2 
This theorem states that a Mellin transform of gives rise to a family of zeta 

functions. In this section we first derive the terminating hypergeometric series for the 
polynomial factor of Eq. (2). We have 
Proposition 1. For q a nonnegative integer 

(-1)1 



P q {s) = H q (s) + 



~ 1)q F ( a S - 1 -2 
-p 2 F^-g,-,-,2 



(5) 



Once this form of P q {s) is obtained, one may determine the functional equation 
P q (s) = {—\) q P q {l — s) and that its zeros are all simple and lie only on Re s = 1/2 

IE]. 

Proof of Proposition 1. We begin by writing 

_l)k 2 3(q-k) ( s 



H q (s) 



9-1 

E 



(6) 



, k\(2q-2k)\ \2j q /, 
and using l/(x - n)\ = (-l) n (-x) n /T(x + 1), {x) n - k = (-l) fc (x) n /(l - x - n) k , and 
the duplication formula (x)2 n = 2 2n (x/2) n ((x + l)/2)„. These steps give 



H q (s) 



from which we may read off 
(s/2) q 2^ 



{s/2) q 2^ ^ (-g) fc (l/2 - q) k 1 1 
(2q)\ ^ (l-s/2-q) k k\2^ 



(7) 



H q {s) 



(2(7)! 



F 1 (-q,l/2-q;l-s/2-q;l/2) 



{-q) q {l/2-q) q 1 1 
(1 - 8 /2 - q) q q\ 2i 



By using (—q) q = (—l) q q\ we then have by re- arranging the last term of this equation 

H q (s) = ^||^ [ 2 Pi(-g, 1/2 -q;l-s/2- q; 1/2) 
l\ g sin(7rs/2) r(s/2) 



-r(l-s/2-g) 



(9) 



r(l/2-g) 

The 2 -Pi function of Eq. (9) can be transformed to another hypergeometric function 
at argument —1/2, —1, or 2. In particular, if we apply ([15], 9.132.2) we have 

2 F 1 (-q, 1/2-q; l-s/2-g; 1/2) = Q|±^r(l-s/2-g) sin(vr S /2)r( S /2) 2 F x {-q, s/2; 1/2; 2) 

(10a) 

(_Jl ' ,,(f/ l/2 ^r( S /2) 2J F\(-g, S /2;l/2;2). (106) 



^T(s/2 + q) 

In obtaining Eq. (10a) we used T(l/2 — g)T(l/2 + q) = n / cos(irq) = (— 1)% for q a 
nonnegative integer. By then inserting Eq. (10b) into Eq. (9) and re-arranging we 
have 

H q (s) = ( 2 g ), r(9 ^ F 1/2) (- 1 ) 9 t ^(-9> »/2; V2; 2) - 1] • (11) 
We then use T(g + 1/2) = ^(2q - l)!!/2 9 and (2q - 1)!! = {2q)\/2 q q\ to arrive at 

H q (s) = [ 2 Pi(-g, s /2; 1/2; 2) - 1] , (12) 

giving Eq. (5). 

Corollary 1. By applying the transformation formula [15] 2 Fi (a,/3; 7; z) = 
(1 — z)~ a 2 F\ (a, 7 — (3; 7; we immediately find the functional equation P q (s) = 
(_l)«P fl (l- s ). 

Corollary 2. From the functional equation we have P^(s) = (— 1) 9+3 °P^(1 — s) 
so that P^')(l/2) = when q + j is an odd integer and in particular P 9 (l/2) = when 
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q is odd. The latter fact may also be found from the a = — 1/2 case of the answer to 
the following question. For what values of n is 2 Fi(—n, (a + l)/2; a + 1; 2) = 0? The 
hypergeometric function here may be written as an nth divided difference summation 
and as a result 2 F x {-n, {a + l)/2; a + 1; 2) = r(l/2)r[-(n + a)/2]/r(-a/2)r[(l - 
n)/2]. The simple poles of the last T factor dictate that the 2 F\ function vanishes 
when n is an odd integer, as expected. 

The following subsumes Theorem 3.4 of Keating [T7j for special values of the 
function C, 2 qfi{s). 

Proposition 2. Put c q = (— l) q (2q)\/q\ = H 2q (0) and let Bj denote the Bernoulli 
numbers. For n, q = 0, 1, 2, . . . and m = 1, 2, ... we have 

(i) C 2g ,o(s) = c, 2 F 1 (-q, -s/2; 1/2; 2)C(s), 

(ii) C 2g ,o(2m) = c, 2J P 1 (-g, -m; 1/2; 2)(2tt 



2(2m)\ 

hi) C2,,o(-n) =c, 2 Fi(-g, n/2; 1/2; 2)- " ' 



-l) n B n+ i 



and 

( iv ) C 2(? ,o( S ) = C 9 

where 



71 + 1 



2 F 1 (-g, s/2; 1/2; 2)C'(s) + C{s)^- ^(-q, s/2; 1/2; 2) 



^ 2 F,(- 5 , s/2; 1/2; 2) = 5 E {t/|+/ 2 >J W'< s / 2 + ■>> " *( s / 2 " I 
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and ip = T'/r is the digamma function [2j [TJ [15]. In particular, we have the values 
((0) = -1/2, C'(0) = -ihi27r and 



v 



^ 2 F 1 (-q,s/2; 1/2; 2) 



^ E TTT^TT = ^ 3^(1 - Q, 1, 1; 3/2, 2; 2). (14) 



s=0 -ft (1/2); j 



Proof of Proposition 2. Parts (i), (ii), and (iii) are obvious from Eqs. (3) and (5) 
and the value of the Riemann zeta function at the negative integers or at the positive 
even integers. Part (iii) follows by writing the series form of the function 2-F1 and using 
the derivative of the Pochhammer symbol (d/dz)(a) n = (a) n [ip(a + n) — i[>(a)]. The 
second line of Eq. (13) follows by applying the functional equation of the digamma 
function [T5] . 

Part (v) can be obtained as the limit s — > of the result of (iv). We provide 
another proof using a representation of the Pochhammer symbol in terms of Stirling 
numbers of the first kind s(j, k) pQ. We have 

n n 

(z) n = Y[(z + k — l) = ^(-l)"- fe S (n, k)z k . (15) 

k=l k=0 

Therefore we have (d/dz)(z) n \ z= o = (n — l)\, where we used s(n, 1) = (— l) n_1 (n — 1)!. 
Using this fact in the series form of 2 -Pi gives the first equality in Eq. (14). For the 
second equality in Eq. (14) we first shift the summation index in the previous one: 



d 2iM-<2,s/2;l/2;2) 



ds 



g<-«W ^ . (16) 

„„ H (V2)i+i U + 1) 



We then use l/(j + 1) = (l)j/(2)j, apply the property (a)j+i = a(a + and the 
rest of Eq. (14) follows. 
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Remarks. Part (iii) of the Proposition covers the trivial zeros of the zeta function 
when n = 2m is an even integer. In regard to part (iv), ('(s) is often easily found in 
terms of ((s) itself from the functional equation. For instance, we have £'(— 2n) = 
(-l) n (2n)\((2n + l)/2(2vr) 2n for n = 1,2,.... By the method of part (iv) higher 
derivatives of (2 q ,o( s ) ma Y be computed. 

Representation of the Riemann xi function 

Similarly to the £ = case of Eq. (1), we put 6j(x) = Y^L-oc f^j{ n V^) where 
f n (x) = (8TT)- n / 2 H n (V2^x)e" wx2 . We then put 

1 oo 

^(x) = -[Oi{x) - / 2 ,(0)] = £ f2j(nVx), (17) 

Z n=l 

where explicitly [15] /^(O) = (— l)- ? (47r)~- ? (2j — 1)!!. The functional equation of 6j 
carries over to ipj so that the following fact is essentially given in Ref. [5]. 
Lemma 1. The function ipj satisfies 









(i) 










\[x 



/2i(0). (18) 

We let as usual £ be the Riemann xi function, defined for all complex s as £(s) = 
s(s — l)ir~ s / 2 T(s/2)((s) /2 and satisfying £(s) = £(1 — s). Very closely related to 
Eq. (5), we put p n (s) = (87r)-"(-l) n (2n)! 2 F l (-n, s/2; 1/2; 2) /n\. We then have 
Propositions 3 and 4. 

Proposition 3. For b > there holds 

Pj {s) = f°° x- {s+1)/2 ^{x)dx + f°° tf> j (x)x a/2 - 1 dx 

S(S — 1) Jl/b Jb 



(19) 



Proposition 4. For complex b with |arg b\ <n/2 we have (i) 



F 6 (s) + (-l)'F 6 -i(l-s) 



(20) 



and (ii) when additionally \b\ = 1 we have 



Pi(s) 



F 6 (s) + (-l)'F 6 (l-s), 



(21) 



where we define 



F b (s)= / ^Wx'^dx-faip)- 

Jb 



(22) 



When j = 0, /o(0) = 1 and Proposition 3 reduces to the original equation of 
Riemann [27], while for the special case b = 1 a result of Bump and Ng [5] and 
Keating [TTJ (Theorem 3.2) is recovered. 

For the proof of Proposition 3 we alternatively evaluate the integral 



as Jo°° ij)j{x)x a l 2 ~ l dx = Jq ip j (x)x S//2 ^ 1 dx + ip j (x)x S//2 ~ 1 dx. We then put y = 1/x in 
the integral on [0,6], apply Lemma 1, and the result Eq. (19) follows. 

Proof of Proposition 4. (i) Equation (19) holds not just for positive b but for all 
values b in the wedge |arg b\ < it/ 2 where the theta function ipj(x) is well defined, 
(ii) The complex conjugate of the function F b (s) of Eq. (18) is F^(s) so that Eq. (21) 
holds whenever b = b~ l . That is, part (ii) holds whenever b lies on the unit circle 
between —i and %. 




(23) 
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Final remarks 

(i) The speciality of the points x = ±i for the function ipj(x) can be seen by 
writing 

oo oo 

$j(x) = (8n)- j H 2j (V2^n)e- nn2x = (8vr)- J ' £ H^V^n^™ 2 ^-™ 2 ^ 

n=l n=l 



(8tt)- j Y, ^2i(v^7ran)(-l) n e-" n M 

n=l 



= -(8ir)- j ]T H 2j (V2^x~n)e-™ ^ + (8n)- J ]T H 2j {V2^2k) e - 4nk M . (24) 

n=\ k=l 
nodd 

(ii) The properties of a theta function have also been useful in an inverse scattering 
approach to the Riemann hypothesis (e.g. [9]). (hi) Proposition 4 offers a prospect 
for generalizing the very important Riemann- Siegel integral formula [12] that takes 
the form 

where F(s) = n- s / 2 T(s/2)I(s) and I(s) is a certain contour integral taken over a 
straight line of slope —1 crossing the real axis between and 1 and directed from 
upper left to lower right. 
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